SPECTRUM OF PERMANENT'S VALUES AND ITS 
EXTREMAL MAGNITUDES IN A 3 n AND A n (a,f3,j) 



g ■ VLADIMIR SHEVELEV 

W)' Abstract. Let A^ denote the class of (0, 1) square matrices containing 

in each row and in each column exactly k Vs. The minimal value of fc, 

■ for which the behavior of the permanent in A^ is not quite studied, is 
k = 3. We give a simple algorithm for calculation upper magnitudes of 
permanent in A^ and consider some extremal problems in a generalized 
class A„(a, f3, 7), the matrices of which contain in each row and in each 

■ column nonzero elements a, /3, 7 and n — 3 zeros. 
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1. Introduction 



Let A„ denote the class of (0,1) square matrices containing in each row 
and in each column exactly k l's. If A 6 A 3 , then matrix k~ x A is doubly 
stochastic. Therefore, A^-matrices are also called doubly stochastic (0,1)- 
matrices (cf. [II]). Furthermore, for a given real or complex numbers 
a, /3, ... 7, denote A n (a, /?,... 7) the class of square matrices containing ev- 
ery number from {a, ...7} exactly one time in each row and in each 
column, such that the other elements are 0's. 

•l-H ' 
X ■ 

Definition 1. We call p-spectrum in A^ (denoting it ps[A^\) the set of all 



the values which are taken by the permanent in A 



Note that p-spectrum in A* trivially is {1}. It is known (cf. Tarakanov 
1251) that 



p S [A^] = {2,2 2 ,2 3 ,...,2LfJ}. 

But, for k > 3, p-spectrum of A^, generally speaking, is unknown. Green- 
stein (cf. [IT] , point 8.4, Problem 3) put the problem of describing the 
p-spectrum in A 3 . In this paper we find p-spectrum on symmetric matrices 
in A 3 with ones on the main diagonal and give an algorithm for calculation 
upper values of p-spectrum in A 3 . We also obtain several results for a gen- 
eralized class A n (ct, f3, 7) with real nonzero numbers a, (3, 7. Some results of 
the present paper were announced by the author in 
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2. What is known about A^? 

1) Explicit formula for \A 3 n \ (cf. Stanley [21], Ch.l) 

(-l) k2 n\ 2 {k 2 + 3k 3 )\2 kl 3 k2 



(2-1) |A£| = 6-» £ 



2) Asymptotic formula for \A 3 n \ (cf. O'Neil [12]) 

(2-2) ^l = ||^e- 2 (l + 0(n-^)), 

where e > is arbitrary small for sufficiently large n. 

In addition, note that |A n (a, (3, j)\ with different a, (3, 7 is, evidently, the 
number of 3-rowed Latin rectangles of length n such that 

(2.3) |A n (a,/3, 7 )| = n\K nj 

where K n is the number of reduced 3-rowed Latin rectangles with the first 
row {l,2,...,n}. It is known (Riordan [13J , pp. 204-210) that 

LSJ , x 

(2.4) K n = ( J D n „ k D k U n - 2k , 
where D n is subfactorial. 

(2.5) D = 1, D n = 71ZV1 + (-l) n , n > 1, 

is sequence of Lucas numbers of the Menage problem which is defined 
by Cayley recursion (cf. [13], p. 201) 

E/b = l, ^i = -l, ^2 = 0, 

n (—l) n 

(2.6) £/ n = nf/„_ 1 + -C/ n _ 2 + 4^ n>3 

n — 2 n — 2 

(see |23l, sequences A102761, A000186). 

Denote, furthermore, A n the set of matrices in A^ with l's on the main 
diagonal. Note that 

(2.7) ps[A 3 n } = ps[Kl]. 

Indeed, it is well known that every A^-matrix A has a diagonal of ones (i.e., 
a set of l's no two in the same row or column). Let I be such a diagonal. 
There exists a permutation of rows and columns n such that will be the 
main diagonal of tt(A). Nevertheless, per(^(A)) = per A and (12. 7p follows. 

3) A known explicit formula for |A n | (Shevelev [19]) has a close structure 
to ([23D: 

I - 1 

3 fn\ 

(2.8) |AJ = ( \S n -kSkU n -'2i: 



fc=0 
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where sequence {S n } is defined by recursion 

(2.9) 50 = 1, 51 = 0, S n = (n- 1) (S n -i+ \s n „ 2 ), n>2. 

4) Asymptotic formula for |AJ (Shevelev [19J) 

(2.10) \-Xl\ = cV^(-) 2n (l + 0(n- 1+£ ), 

e 

where 

C = 2Vire- 5 = 0.29098... 

and e > is arbitrary small for sufficiently large n. 

5) Denote the set of symmetric matrices in A^. P-spectrum on A^ is 
given by the following theorem (Shevelev [16J) 

Theorem 1. Let R(n;3) denote the set of all partitions of n with parts 

more than or equal to 3. To every partition r e R(n; 3) : n = n\ + n 2 + ... + 

n m , m = m(r), put in a correspondence the number 

in 

(2.11) H(r) = l[a(n i ), 

i=l 

where sequence {a(n)} is defined by the recursion 

(2.12) a(3) = 6, a(4) = 9, a(n) = a(n - 1) + a(n - 2) - 2, n > 5. 
Then we have 

(2.13) ps[Al] = {H(r):reR(n;3)}. 

6) The maximal value M(n) of permanent in A^ was found by Merriell 

Theorem 2. If n = h (mod 3), h = 0, 1,2, then 

(2.14) Af(n) = 6 a S S L(|) fc J- 

Note that, the case h = of (12.141) easily follows from a general Minc- 
Bregman inequality for permanent of (0,l)-matrices (see [IT] , point 6.2, and 
[4])- 

7) Put M(n) = M^){n). In case of n = (mod 3) , Bolshakov [3] showed 
that the second maximal < M^)(n) of permanent in A^ (such that 
interval (M^ , MW) is free from values of permanent in A^) equals to 

f 20, i/n = 6, 

(2.15) M (2) (n) = < 120, if n = 9 

|j6f, i/ n>12. 

Note that both M (1) (n) and M (2) (n) are attained in A^ (Shevelev [IB]). 
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8) Denote m(n) the minimal value of permanent in A 3 . In 1979, Voorhoeve 
[26] obtain a beautiful lower estimate for m(n) : 

(2.16) m(n) > 6(-) n " 3 . 

3 

This estimate remains the best even after proof by Egorychev [7] and Falik- 
man [S] the famous Van der Waerden conjectural lower estimate per A > 
for every nxn doubly stochastic matrix A. Indeed, this estimate yields only 
m(n) > 3™^-, such that (12. 161) is stronger for n > 4. 

9) Bolshakov [2] found p-spectrum in A 3 in cases n < 8. Namely, he added 
to the evident p-spectrums ps[A|] = {6} and ps[A|] = {9} also the following 
p-spectrums 

ps[Al] = {12, 13}, ps[Al] = {17, 18, 20, 36}, 

(2.17) ps{A*} = {24, 25, 26, 27, 30, 31, 32, 54}, 

ps[Al) = {33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 44, 45, 48, 49, 52, 72, 78, 81}. 

3. A GENERALIZATION OF THEOREM 1 ON MATRICES OF CLASS 
A n (cu, /3, 7) WITH SYMMETRIC POSITIONS OF ELEMENTS 

Denote A n (a, j3, 7) the set of matrices in A n (a, (3, 7) with /3's on the main 
diagonal. It is clear that, together with (12.31) . 

(3.1) \K{a,/3,j)\ =K n . 
Note that, as for sets A 3 , A n , we have 

(3.2) ps[A n (a, f3, 7)] = ps[A n (a, (3, 7)]. 

Denote, furthermore, A n (a, (3, 7) the set of matrices M = {rriij} in A n (a, f3, 7) 
with symmetric positions of elements: m^j = a if and only if rrijj = 7. 
P-spectrum on A n (a, /?, 7) is given by the following theorem. 

Theorem 3. // to every partition r G R(n; 3) : n = rii + n 2 + ... + n m , m = 

m(r), corresponds the number 

(3.3) H a>/3>J (r) = JJa(rij), 

8=1 

where sequence {a(n) = a(a,/3,j; n} is defined by the recursion 

a(3) = a 3 + /3 3 + 7 3 + 3a^7, 
a(4) = a 4 + /3 4 + 7 4 + 4a/3 2 7 + 2(a 7 ) 2 , 
a(n) = f3a(n — 1) + crya(n — 2) + 

(3.4) a"" 1 (a- / 5-7) + 7 "- 1 (7-/3-a), n > 5, 



SPECTRUM OF PERMANENT'S VALUES IN A^ AND A n (a,/3,j) 



5 



then we have 



(3.5) 



ps[A n (a,(3,j)} = {i/ ai/3j7 (r) : r G i?(n;3)}. 



Proof. Let S n be the symmetric permutation group of elements 1, ...,n. 
Two positions ji), (12,32) are called independent if ^ 7^ jfe, k=l,2. We 
shall say that in the n x n matrix M = {m^} a weight is appropriated 
to the position Let s £ S n has not any cycle of length less than n. 

Consider a map 



appropriating to the position s J (l) the weight m^. 

Lemma 1. 1) the map o is bijective; 2) if E is a set of pairwise independent 
positions, then cr(E) is also a set of pairwise independent positions. 

Proof, a) Consider two distinct positions 



Let i\ 7^ Z2 such that, say i\ > i 2 . Show that s n (l) 7^ s* 2 (l). Indeed, if to 
suppose that s ll (l) = s l2 (l), then s ll ~ 12 = 1, i.e., s has a cycle of length 
h — 22 < n in spite of the condition. Conversely, if s n (l) 7^ s* 2 (l), then 
&i 7^ since s -1 has not any cycle of length less than n as well, 
b) Let positions (13. 6p be independent. The both of inequalities (13. 7p hold 
and, as in a), we have s n (l) 7^ s* 2 (l), s jl (l) 7^ s j2 (l), i.e. the positions 
&((h, Ji ) ) , ji)) are independent as well. 
■ 

Lemma 2. Lei s G S* n /iai>e not any cycle of length less than n. Then 
(0, l)-matrix S having l's on only positions 



is a incidence matrix of s. 

Proof. Since s has not cycles of length less than n, then {s 1 (l), s n (l)} 
is a permutation of numbers {l,...,n}. Thus the set of positions of l's 
of matrix S coincides with the set of l's of the incidence matrix of s : 



Let P = P n be (0, l)-matrix with l's on positions (1, 2), (2, 3), (n — 
l,n), (n, 1) only. 




(^(1), , 2 (1)), (, 2 (1), , 3 (1)), (^(l), ,"(1)), (s n (l), s\l)) 



(l,s{l)),...,{n,s(n)). 
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Lemma 3. Let s G S n have not any cycle of length less than n. If S and 
S^ 1 are the incidence matrices of s and then we have 

(3.8) <r-\S) = P, o-- 1 ^- 1 ) = P- 1 . 
Proof. Both of formulas follows from Lemma [2j 

■ 

Noting that a (I) = I, where / is the identity matrix, we conclude that 

(3.9) S- 1 + I + S = a{p- 1 + I + P). 

Moreover, since, by the bijection a, to every diagonal (i.e., to every set of n 
pairwise independent positions) of the matrix aS^ 1 + (31 + ^yS corresponds 
one and only one diagonal of the matrix ctP -1 + (31 + 7P with the same 
products of weights, then we have 

(3.10) per(aS~ 1 + (31 + 7S) = per(aP- 1 + (31 + jP). 

Note that from the definition it follows that, for every matrix M G A n (a, (3, 7), 
we have a representation 

(3.11) M = aS' 1 + (3I n + jS, 

where S is the incidence matrix of a substitution s. In case when s has not 
any cycle of length less than n, the matrix M is completely indecomposable 
matrix in A n (a, (3, 7). Thus, by (13.101) . all completely indecomposable matri- 
ces of A n (a, (3, 7) have the same permanent, equals to per(aP~ 1 +/3I n + 'yP). 

In general, a substitution s with the incidence matrix S in (13. lip can- 
not have cycles of length less than 3. Indeed, if for some i, we have either 
s(i) = i or s(s(i)) = i, then in both cases s(i) = s~ 1 (z) which means coinci- 
dence of positions l's of the matrices S and S~ l in the z-th row. 

Let s G S n be an arbitrary substitution with cycles of length more than 
2. Let 

r 

(3.12) s = Y[ Sj , 

3=1 

where Sj G S^, lj > 3, Y7j=i h = n -> be the decomposition of s in a product 
of cycles. Then the matrix M = aS" 1 + (3I n + 7^ is a direct sum of the 
matrices Mj = cuS^ 1 + ^I + 'yS^ such that, by (I3.10p . perMj = per(aP^ 1 + 
(3Ii 3 + 7P) and we have 

r r 

(3.13) perM = JJperM,- = JJper(aP _1 + (3I h + 7P). 

3=1 3=1 

It is left to notice that Mine [10] found a recursion (13.41) for per(al n + (3P + 
7P 2 ) and, as well known, the multiplication an n x n matrix by P _1 does 
not change its permanent. 
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Therefore, per(aP~ l + fH\. + jP) = per(al n + j3P + jP 2 )- 

■ 

Example 1. Let us find ps[A n (— 1, 3, 2)]. 

We have the following partitions of 11 with the parts not less than 3 

11 = 8 + 3 = 7 + 4 = 6 + 5 = 3 + 4 + 4 = 3 + 3 + 5. 

According to (13. 4p . for a(n) = a(— 1, 3, 2; n), we have a (3) = 16, a (4) = 
and for n > 5, 

a(n) = 3a(n - 1) - 2a{n - 2) + 6(-l) n . 
Using induction, we find 

i/ n is odd, 



am 



2 n+1 + 2, i/ n is even. 



Therefore, 

pstAnC-1,3,2)] 



{a(ll), a(3)a(8), a(5)a(6), a 2 (3)a(5), a(3)a 2 (4)} = 
{4096, 8224, 8320, 8704, 16384, 18496}. 

■ 

In the following examples we calculate p-spectrum for arbitrary n. 
Example 2. Let us find ps[A n (— 1, 2, 1)]. 

By induction, for a(n) = a(— 1, 2, 1; ra), we have 

I 2, if n is odd, 

a \ n ) = \ A - f 

14, z/ n zs even. 

Further, again using induction, one can find that, if n is even, then 

ps[A n (-l,2,l)] = {4, 4 2 ,...,4L(f>J} 

and, if n is odd, then 

ps[A n (-l,2,l)] = {2, 2-4, 2-4 2 ,...,2-4^)J}. 



Example 3. Analogously, in case o/A n (— 1, 1, 1), fora(n) = a(— 1, 1, 1; 
we /iave 
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4, if n = (mod 6), 
a(n) = ^ —2, if n = 3 (mod 6), 
1, otherwise. 



and 



p S [A„(-l,l,l)] = 
{1, -2, 4,...,(-2)L(^)J}, i/ n = 1,2 (mod 3), 
{1, -2, 4,..,(-2)L(t)J ) (— 2)L(t)J }, if n = (mod 3). 
It is interesting that, in case of n multiple of 3, the permanent omits the 
value (-2)L(^)J. 
■ 

4. Merriell type theorems in a subclasses of A n (a !( 5,7) 

Note that in class A n (a, f3, 7) the Minc-Bregman inequality and the Mer- 
riell theorem , generally speaking, do not hold even for positive a, (3, 7. 
Nevertheles, some restrictions on a, f3, 7 allow to prove some analogs of 
the Merriell theorem. Recall that M(n) ( 12.141) is attained in A 3 . Denote 
M n (a, f3, 7) the maximal value of permanent in A n (a, f3, 7). 

Theorem 4. Consider a class A n (a,/3, 7) with the numbers a, 0, 7 satisfy- 
ing "triangle inequlities" 

(4.1) 0<a<,5 + 7, 0<7<a + /3, 
and the following additional conditions 

(4.2) a 3 (4) < a 4 (3), a 7 + /3(a(3))s < (a(3))§, 

where sequence {a(n)} zs defined by recursion < \3.4[ - Then, forn multiple of 
3, we have 

(4.3) M„(a,/3, 7 ) = (a(3))i. 

Proof. Note that conditions (I4.ip - fl4.2p are satisfied, e.g., in case a = 
(3 = 7 = 1. Using induction, let us prove that 

(4.4) a(n) < (a(3))t. 

Indeed, for n = 3, this inequality is trivial, while, for n = 4, it follows from 
the first condition f)4.2p . Let it hold for n < m — 1. Then, according to (13. 4p . 
we have 

a(m) = pa(m - 1) + a 7 a(m - 2) + a m ~\a - - 7) + 7 m " 1 (7 - at - 0) < 

/3(a(3))— +a 7 (a(3))— = 
(a(3))^(/?(a(3))3 + 07) < (a(3))^(a(3))i = (a(3))f . 
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Note that, according Theorem (3J the equality in (14.4p holds in a direct 
sum of (3 x 3)-matrices of A 3 (a,/3,7) which corresponds to the partition 
n = 3 + 3 + ... + 3. Let now A G A„(a, /3, 7). By Theorem [31 there exists 
a partition of n with the parts not less than 3, n = n± + ... + n m , such that 

rn 

per A = Y\ a>(n>i) 

i=3 

and, in view of ( 14. 4ft . we have 



per A < JJa(3)^ = (a(3))f. 



8=3 



This proves ( 14. 3ft . 



Example 4. Consider case (3 = 7 — a. 

Let us find the values of a, depending on the magnitude of 7, for which 
the conditions of Theorem H] are satisfied. According to (13. 4p , we have 

(4.5) o(3) = a 3 + (7 - a) 3 + 7 s + 3«( 7 - 0)7 = 2 7 3 , 

(4.6) a(4) = a 4 + (7 - a) 4 + 7 4 + 4a(7 - a) 2 7 + 2(aj) 2 = 2(a 4 + 7 4 ). 
Thus the condition a 3 (4) < a 4 (3) means that a 4 + 7 4 < 237 4 , or 

(4.7) < a < (25 - 1)*7 = 0.7140199. ..7. 

and it is easy to verify that the second condition in (14. 2 j) is satisfied as well. 
As a collorary, we obtain the following result. 

Theorem 5. 7/ ( |^ . 7 ) holds, then, forn multiple of 3, we have 



(4.8) M n (a, 7 -a,7) = 2t 7 ". 

■ 

Simple forms of sequence {a(n)} in Examples [UTS] allow to suppose that 
in case (3 = 7 — a (or symmetrical case /3 = a — 7) sequence {a(n)} keeps 
a sufficiently simple form. We find this form in the following lemma. 

Lemma 4. If [3 — 7 — a, i/ien sequence {a(n)} which is defined by recursion 
<\3-4\) has the form 

2j n , if n is odd, 

2(a n + 7 n ), if n is even. 



(4.9) a(n) 



Proof. Using induction with the base (14. 5 ft - (14.61) . suppose that (14.9ft 
holds for m < n. Then, by H3.4[) . for even n, we have 

a{n + 1) = (7 — a)a{n) + a^a{n — 1) + 2a n (a — 7) = 



SPECTRUM OF PERMANENT'S VALUES IN A^ AND A n (a,/3,j) 



10 



2( 7 - a){a n + 7") + 2a 7 n + 2a n {a - 7) = 2 7 n+1 , 
while, if n is odd, then we have 

a(n + 1) = 2(7- a) 7 n + 
2a 7 (a"- 1 + 7"" 1 ) + 2a"(a - 7) = 2(a" +1 + 1 n+l ). 

■ 

Let now 

(4.10) a > (23 - 1)27 = 0.7140199. ..7. 

Theorem 6. // ( |^ . i (J\) holds, then, for n multiple of 4, we have 



(4.11) M„(«, 7 -«, 7 ) = (2(« 4 + 7 4 )) f . 
Proof. From p~5]) , fT4TB]) and (14.101) we conclude that 

(4.12) a(3) < (a(4))i. 
Let us show that, for n > 3, 

(4.13) a(n) < (a(4))l 

For n = 4, inequality (14. 13j) is trivial. For n > 5, we have 

a" + 7" = a n (l + (^D < a n (l + (^) 4 )? < (a 4 + 7 4 )^ 
and thus, using Lemma HI we have 

a(n) < 2(a n + 7 n ) < 2%(a 4 + 7 4 )f = (a(4))?, n > 3. 

Let now A G A n (a,/3, 7). By Theorem [3j there exists a partition of n with 

the parts not less than 3, n — n\ + ... + n m , such that 

m 

per A = Y\ a(ni) 

i=3 

and, in view of ( I4.13p . we have 

m 

nil: n 
a(4)- = (a(3))4 

i=3 

with the equality in a direct sum of (4 x 4)-matrices of A 3 (a,/3,7) which 

corresponds to the partition ?t, = 4 + 4+ ...+4. 

■ 

Note that, if a 7^ (23 — 1)27, then in Theorem 5 we have only maximizing 
matrix (up to a permutation of the rows and columns) which corresponds to 
the partition n = 3 + 3 + ... + 3;in Theorem [6] we also have only maximizing 
matrix (up to a permutation of the rows and columns) which corresponds 
to the partition n = 4 + 4 + ... + 4. It is interesting that, only in case of the 
equality 
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a = #7, where 6 = (2a — 1)?, the both of Theorems 5-|6]are true for every 
n multiple of 12 with the equality of the maximums: (2(ce 4 + 7 4 ))^ = 2 §7™. 
Thus, up to a positive factor 7, tie ciass 

(4.14) An(M-M), = (2? -1)2, 

possesses an interesting extremal property: it contains f| + 1 maximizing 
matrices (up to a permutation of the rows and columns), instead of only 
maximizing matrix, if 9 ^ (2 s — 1)4. 

Indeed, the number of the maximizing matrices (up to a permutation 
rows and columns), is defined by the number of the following partitions of 
n = (mod 12) : 

n = 3 + 3 + ... + 3, n = (4 + 4 + 4) + (3 + ... + 3), 

n 

n = 4 + ...+4 + 3 + ... + 3, i = 0, 1, 

3i n—VM, 

■ 

5. Estimate of cardinality of p-spectrum on circulants in 

and A n (a,/3, 7 ) 

Denote C the set of the circulants in A^. Note that a circulant 
A e Al has a form A = P l + P j + P k , 0<i<j<k<n, where P = P n 
is (0, l)-matrix with l's on positions (1, 2), (2, 3), (n — l,n), (n, 1) only. 
Multiplicating A by P~ l we obtain circulant B of the form 

B = I n + P r + P s 

with peri? = per A. Since I? is defined by a choice of two different values 
< r < k < n, then trivially 

< (;) < f . 

Now we prove essentially more exact and practically unimprovable estimate. 
Theorem 7. VKe /iai>e 

(5.1) Ps[K)<l^\- 

Proof. Let us back to the general form 

A = P { + P j + P k , < % < j < k < n. 

Note that, A is defined by a choice of a vector k), but its rotation, i.e., 
a passage to a vector of the form (i + l, j + 1, k + l) (mod n) , does not change 
the magnitude of per A. Indeed it corresponds to the multiplication A by 
P , and our statement follows from the equality per(P l A) = per A. Besides, 
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its reflection relatively some diameter of the imaginary circumference of 
the rotation, by the symmetry, keeps magnitude of the permanent. Since 
geometrically three points on the imaginary circumference define a triangle, 
then our problem reduces to a triangle case of the following general prob- 
lem, posed by Professor Richard H. Reis (South-East University of Mas- 
sachusetts, USA) in a private communication to Hansraj Gupta in 1978): 

" Let a circumference is split by the same n parts. It is required to find the 
number R(n, k) of the incongruent convex &-gons, which could be obtaind 
by connection of some k from n dividing points. Two /c-gons are considered 
congruent if they are coincided at the rotation of one relatively other along 
the circumference and (or) by reflection of one of the /c-gons relatively some 
diameter." 

In 1979, Gupta [6] gave a solution of the Reis problem in the form (for a 
short solution, see author's paper [22J): 



Z V LUJJ / K d|g cd (fc,n) X ~ d 

If to denote A^ C Aj; the set of the circulants in A^, then from our arguments 
it follows that 

(5.3) ps[A k n }<R(n : k). 

In case k = 3, from fl5.2l) - fl5.3p we find 

' x§j if n = (mod 6), 
— jyi, if n = 1,5 (mod 6), 
if n = 2,4 (mod 6), 
™ 2+3 if n = 3 (mod 6), 



p4K\ < { 

and (15. ip follows. 



12 ' 



Example 5. In case n = 5 we have only two incongruent triangles corre- 
sponding to circulants 1$ + P + P 2 and 1$ + P + P 3 . 

Nevertheless, the calculations give per (I 5 + P + P 2 ) = per (I 5 + P + P 3 ) = 
13. Thus ps[Al] = {13}, and bs[A 3 ]| = 1. 

Example 6. In case n = 6 we have three incongruent triangles correspond- 
ing to circulants I 6 + P + P 2 , I e + P + P 3 and I 6 + P 2 + P 4 . 

The calculations give per(I 6 + P + P 2 ) = 20, per(I 6 + P + P 3 ) = 17, 
while per(I 6 + P 2 + P 4 ) = 36. Thus ps[A 3 6 ] = {17, 20, 36}, and |ps[A 3 ]| = 3. 
Note that a respectively large magnitude of per{I & + P 2 + P 4 ) is explained 
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by its decomposability in a direct product of circulants (J3 + P + P 2 ) <g> 
(J3 + P + F 2 ), such that per(J 6 + F 2 + P 4 ) = (per(I 3 + P + P 2 )) 2 = 6 2 = 36. 

It is clear that, in case of circulants in A n (a, (3,j), the upper estimate 
( 15. II) yields either the same estimate, if a = j3 = 7, or L^ij^J , if ct = /? 7^ 7 
(and in the symmetric cases), or [^^J > if a i/^7 are distinct numbers. 

Add that a bijection indicated in [22] allows to apply formula ( 15. 2 j) to 
enumerating the two-color bracelets of n beads, of which are black and n — 
k are white (see, e.g., the author's explicit formulas for sequences A032279- 
A032282, A005513-A005516 in [23J). 

6. Algorithm of calculations of upper magnitudes of the 

permanent in a^ 

Theorem [Hallows, using some additional arguments, to give an algorithm 
of calculations of upper magnitudes of the permanent in A^. In connection 
with this, we need an important lemma for numbers (I2.12p . 

Lemma 5. For ni,n 2 > 3, we have 

(6.1) a(rii + n 2 ) < a(ni)a(n 2 ) ,rzi,n 2 >3. 

Proof. By usual way, from (I2.12p we find 

(6.2) a(n) = <p n + 2 + {-l) n tp- n , 
where (p = ^ is the golden ratio. 

Denote e(n) = (—l) n ip~ n . Since n > 3, then \e(n)\ < 0.24, and, conse- 
quently, if n — n\ + n 2l then (2 + e{n\))(2 + e{n 2 ) > 1.76 2 > 3. Therefore, 
we have 

0(7*1)0(712) = (^ ni + 2 + < n i))(v n2 + 2 + e{n 2 )) > 

(6.3) ip ni+n2 + 3 > cp ni+n2 + 2 + e(m + n 2 ) = a(m + n 2 ). 
■ 

Note that, actually, the difference between the hand sides of (6.3) more 
than 1.76((^ ni + y? n2 ). 

Let now n = j (mod 3), j = 0,1,2, and t G N. Let R(m\v) denote 
the set of all partitions of n with parts more than or equal to v. For us an 
important role play cases v = 3,4. To r G R(m;3), p G R{m;A) put in a 
correspondence the sets 

(6.4) Hmjzir) = {n rier o r J; H m . 4 (p) = {U Piep a Pi }. 

In case m — 3, when p — 0, let us agree that is a singleton {6}. 
Consider now the set = L\ (n) : 
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4t+j 

(6.5) L? = |J {6^y : y G 11,,. ,,,{,»■ V > 9^6^}. 

Theorem 8. (algorithm of calculation of upper magnitudes of the perma- 
nent in A 3 ). If n > A(3t + j), then the ordered over decrease set gives 
the \l[ \ upper magnitudes of the permanent in A 3 . 

Proof. Note that the proof is the same for every value of j. Therefore, 
let us consider, say, j = 0. From (6.1) it follows that, if r G R(n; 3) contains 
A3 parts 3 and A3 < | — At, then, for y G iJ„_3A 3; 4(p), we have 

6 X3 y < 6f- 4 *9 3t . 

This means that for the formation the list of all upper magnitudes of the 
permanent in A 3 in the condition n > 12t, which are bounded from below 
by 6t _4 *9 3 *, it is sufficient to consider only a part of the spectrum containing 
numbers {6 X:i y}, where y G H n ^,x r ^{p) with the opposite condition A3 > 
I — At . From the equality 3A3 + . . . + n\ n = n with the condition A3 > | — At, 
we have 

4A 4 + ... + n\ n < n - 3(- - At) = 121 

Since 12t does not depend on n, there is only a finite assembly of such 
partition for arbitrary n. This ensures a possibility of the realization of the 
algorithm. 

For the considered r G R(n; 3), for A 3 > | — At, we have H n -,(r) = 6^~~ , 
where y G H m .^(p), and m has the form m = 3i, 1 < i < At. Thus we 
should choose only H n . 5 (r) > 6t _4t 9 3t , and this yields 

?/ > 9 3 *6?- 4 * = g 3 ^" 4 *, i = 1, 2, 4t. 

■ 

In order to use Theorem 8 for calculation the upper magnitudes M^)(n) > 
M^ 2 ){n) > ... of the permanent in A 3 , in case, say, n = (mod 3), 

1) we write a list of partition of numbers 3i, i = 2,3, ...,At with the parts 
not less than 4. 

2) The corresponding values of y we compare with 9 3 *6 i_4 * and keep only 
y > 9 3t 6 i ~ 4t . 

3) After that we regulate over decrease numbers {?/6t _i }. 

Below we give the first 10 upper magnitudes > > ... > M^ 10 \ 
of the permanent in A 3 for n > 24, via numbers {a(n)} (I2.12p . 

{a(3)t = 6t, if n = (mod 3), 
a(4)a(3)^ = §6^, if n = 1 (mod 3), 
a(4) 2 a(3)^ = |6^, if n = 2 (mod 3). 
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Formula ( \Q.Qh shows that M^(n) is attained in A 3 . 

a(4) 3 a(3) : ^ 12 = , i/ n = (mod 3), 

(6.7) = { 0(7)0(3)^ = §6^, if n = 1 (mod 3), 
0(5)0(3)^ = ^6^, i/ n = 2 (mod 3); 

a(6)a(3) £ ar = -6 s , if n = (mod 3), 

(6.8) M (3) = <J o(4) 4 o(3) 2 ^ 1S = §6^, z/ ti ee 1 (mod 3), 



o(8)o(3)tt = §6^, if n = 2 (mod 3); 



n-9 



0(4)0(5)0(3)^ = f6t, if 71 = (mod 3), 

(6.9) M (4) = <j 0(4)0(6)0(3)°^ = §6^, if n = l (mod 3), 
0(4)0(7)0(3)^ = §6^, i/ 71 = 2 (mod 3); 

0(9)0(3)^ = |§6i, n = (mod 3), 

(6.10) M (5) = <J o(4) 2 a(5)a(3)^ = ±§6^, i/ n = 1 (mod 3), 

n — 20 01 _ n — 2 



a(4) 5 a(3)— = §6—, if n = 2 (mod 3); 



o(4)o(8)o(3)^ = ^6f , i/ n = (mod 3), 

(6.11) M (6) = <j a(5) 2 )a(3)^ = |||6t, i/ n = 1 (mod 3), 
a(4) 2 a(6)a(3)^ = §6^, 2/ n = 2 (mod 3); 

a(4) 2 a(7)a(3)^ = §6f , if n = (mod 3), 

(6.12) M (7) = <J o(10)a(3)^ = §§6^, i/ n = 1 (mod 3), 
a(4) 3 a(5)a(3)^ = §6^, if n = 2 (mod 3); 

a(4) 6 a(3) Ii ^ = |^6 5 , n = (mod 3), 

(6.13) = { o(4)o(9)o(3)^ = ||6t , i/ n = 1 (mod 3), 
0(5)0(6)0(3)^ = §§6^, if n = 2 (mod 3); 

n-18 



o(4) 3 a(6)a(3)^ = ^63 , ifn = (mod 3), 

71 — 13 flQV „ n — 1 



) 3 o(6)o(3)^ = A 
(6.14) = { 0(5)0(8)0(3)^ = ^6V, i/ n = 1 (mod 3), 

a(5) 2 a(4)a(3)^ = lg 6^, i/ n = 2 (mod 3); 



0(5)0(7)0(3)^ = ^|6f , 2/ 72 = (mod 3), 

n-19 



403 
1296 

(6.15) M (10) = { a(4) 3 o(7)a(3)^ = §6^, j/ n ee 1 (mod 3), 
o(ll)o(3) 2 T 1 = §6^, if n = 2 (mod 3). 

7. Main conjectural inequality for maximum of permanent in 
completely indecomposable a 3 -matrices 

Denote A 3 1 the set of completely indecomposable matrices in A 3 , i. e., 
the set of A 3 - matrices containing no A^-submatrices. Let \i\ (n) denote the 



SPECTRUM OF PERMANENT'S VALUES IN A^ AND A n (a,/3,j) 16 

maximum of permanent in 1 . Our very plausible conjecture which we 
call "main conjectural inequality (MCI)" is the following. 

Conjecture 1. (Cf. [21] : pp. 165-166) For ri\,n<i > 3, we have 
(7.1) /ui(ni + ra 2 ) < ^i(n 1 )^i(n 2 ). 

In Lemma Owe essentially proved that in subclass A^ the MCI is valid. 

Besides, in all known cases MCI holds. Moreover, as we shall see, our al- 
gorithm of calculation the consecutive upper magnitudes (M = Mi > M 2 > 
...) of permanent in A^, which is based on MCI, reproduces all Merriell's and 
Bolshakov's results for M\ and M 2 . Note also that, for sufficiently large 
n, the number of consecutive upper magnitudes of permanent in A^ grows 
very quickly with every step of extension of the list of known p-specrums for 
small n. E.g., using the found by Bolshakov ps[A?], i < 8, we obtain, for 
sufficiently large n, 4, 7 and 11 upper values of ps[A^] in cases n = 3k, 2k+l 
and 3k + 2 correspondingly. After calculation ps[A<j], the number of upper 
values of, e.g., ps[Ag fc ] increases more than thrice. 

8. Algorithm of calculations of upper magnitudes of the 
permanent in a^ based on mci 

Let n = j (mod 3), j = 0,1,2, and t G N. Let R{m;v) denote the 
set of all partitions of n with parts more than or equal to v. For us an 
important role play cases v = 3,4. To r G R(m;3), p G R{m;4) put in a 
correspondence sets 

(8.1) n m . 3 (r) = {U ner x ri }; 7r m;4 (p) = {U Piep x Pi } 

where x s runs through all values of permanent in set 1 of completely 
indecomposable matrices in A^ ( in case m = 3, when p = 0, let us agree 
that 7r 3; 4 is a singleton {6}). 

Consider now the set E^ = E^(n) : 

u+j 

(8.2) E® = |J {(f * ~y : y G 7T 3i+i;4 (p), y > 9 3 ^6^}. 

i=l 

Theorem 9. (algorithm of calculation of upper magnitudes of the perma- 
nent in A 3 ). If n > 4(3£ + j), then the ordered over decrease set gives 
the \E { t j) \ upper magnitudes of the permanent in A„. 

Proof. We need three lemmas. 



Lemma 6. For n > 4, we have 

(8.3) Pi{n)<3%. 



SPECTRUM OF PERMANENT'S VALUES IN A^ AND A n (a, j3,j) 17 

Proof. Let, firstly, n = (mod 4). Note that /ii(4) = D 4 = 9. Using 
(7.1), we find 

Mi(n) < /ii(4) Ml (n - 4) < ... < (^i(4))? = 3*. 

Let, furthermore, tt. = i (mod 4), z = 1,2,3. Note that, by (I2.17p . ,ui(5) < 
13 < 3 2 ' 5 . Therefore, using (7.1), we have 

^i(n) < (M4)) — (// x (5)) 1 < 3— 3 2 - 51 = 3s. 



Lemma 7. Let n = 3A3+4A4 + ... + nA n be a partition of n with the parts not 
less than 3. If A3 < I, and n has the form n = 31 + 4m, i7ien ; for completely 
indecomposable matrices Ai G Af, z = 3,4, ...,n, we have 

n 

(8.4) JjGoerA) A * < 6 l 9 m . 

i=3 

Proof. Using Lemma El we have 

n 

4A4+...+nA n 1 / — n— A3 / — 4m 

i=3 



Lemma 8. Lei n = 3 A3 + 4A4 + ... + n\ n and 

A 3 < — g-^-4t, n>4(3i + i), 

where t is a nonnegative integer and j is the residue of n modulo 3, j = 
0, 1, 2, i/ien, /or completely indecomposable matrices Ai G Af , z = 3, 4, n, 
we have 



.5) ]J(perA) A * ^ 6^- 4 '9 3 ' +J . 



n—Aj 

perAi)"' < 6 

i=3 

Proof. Put / = ^ - 4t, m = ^ = 3t + j. Now the lemma follows 



3 4 

from Lemma [71 



It is left to note that, after these lemmas, the proof of Theorem [9] is the 

same as proof of Theorem [HI 

■ 

Note that the using of this algorithm is based on the small elements of 
p-spectrum. 

Consider, e.g., case t = 0, j = 2. According to (18.21) . we have 

2 
8=1 
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(8.6) {6^rperA, A e A 3 5 : QperA > 81} U {81 • 6^}. 

Note that, the second set is a simpleton, since, by MCI, /ii(8, 3) < (/ii(4, 3)) 2 
81. Since, by ([2Tfj) . M(5) = 13 < f , then the first set in @£BJ is empty. 
Thus i?Q 2 ^ = Eq{u) is simpleton: 

£f = {81-6^} 

and we have 

M (1) (n) = 81 • 6^, n > 8, 

which corresponds to Merriell's result in case n = 2 (mod 3). 

Further research of the set ( 18. 2ft . using ( 12. 17ft . gives the following results: 

1) j = 0, n > 24. 

M«(n)=6?, M( 2 )(n) = -6f, 

16 

(8.7) M®{n) = ^tf, M< 4 >(n) = ^6*. 

The continuation of this list requires the knowing of ps[Ag]. Note that a 
more detailed analysis shows that after calculation ps[Ag] in this case one 
can obtain the first 12 + \G\ upper magnitudes of the permanent in A 3 , 
where G = ps[A 3 9 ] n ([69, 116] \ {72, 78, 102, 108}). 

2) j = 1, n > 28. 

M«(n) = -6^, M^(n) = -6^, 
2 9 

M ( 3 )(n) = -6^, MW(n) = -6^, 
36 32 

(8.8) M^\n) = -6^, Af«(n) = -6^, M«\n) = — 6^ 

6 15 216 

It is interesting that in this case ps[Ag] is not used up to M^ 7 \ but the 
continuation of this list requires the knowing of ps[A 3 ]. 

3) j = 2, n> 32. 





= j6 . , 




13 n-2 

= — 6 3 

6 


M (3) (n) : 


n-2 

= 2- 6-5-, 


M (4) (n) 


13 n-2 

= — 6~ 

9 


M (5) (n) 


49 n-2 

= — 6~ 
36 




4 n-2 


M (7) (n) : 


31 n-2 

= — 6~ 

24 ' 


M (8) (n) 


81 n-2 

= — 6~ 

64 
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(8.9) M^\n) = V^, M {w \n) = ^6^, M^\n) = 

4 9 32 

Note that the method not only gives a possibility to calculate the upper 
magnitudes of the permanent in A^, but also indicates those direct products 
on which they are attained. E.g., in ( 18. 9 p M 9 is attained on direct products 
of some matrices A; 6 K\ : 

Ah® As® ... ® A,: Ad® A 7 ® A s ® ... ® Ail A® A 4 ® A fi ® A 3 ® ... ® A s . 



n-8 
3 



Note also that the comparison of ( I8.7p -( l8^9|) with (I6.6p -( 16T5|) shows that 
the following calculated are attained in A^, n > 32 : 
in case n = mod 3, 

M (l) 7M (2) iM (3) )M (4). 

in case n = 1 mod 3, 

M (1) 7M (3) )M (4) )M (5) )M (6) )M (7) 

(and is not attained M^ 2 ^); 
in case = 2 mod 3, 

M (l) ; M (2) ^ M (5) j M (7) ^ M (8) ( M (9) j M (ll) 

(and are not attained M^, M^, M^). 

9. Algorithm of a testing the parity of values of the 

permanent in a^ 

It seems that, among all known methods of calculation of the perma- 
nent, only Ryser's method (cf. [11] , Ch.7) could be used for a creating 
an algorithm of a testing the parity of values of the permanent. Let A be 
n x n-matrix. Let A r be a matrix which is obtained by changing some r 
columns of A by zero columns. Denote S(A r ) the product of row sums of 
A r . Then, by Ryser's formula, we have 

perA = J2S(A )-J2S(A 1 ) + 

(9.1) ^5(A 2 )-... + (-ir 1 ^,S(A ri _ 1 )- 

Let now A have integer elements. Introduce the following matrix function 

1, if all row sums of A are odd, 
0, otherwise. 



(9.2) T(A) : 

From (19. ip we have 



perA^^T(A )-^T(A 1 ) + 
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(9.3) £T(A 2 ) + ... + £T(4»-i) (mod 2). 

Using (19.31) . let us create an algorithm of a search of the odd values of the 
permanent in A^. Since, evidently, per A = detA (mod 2), then A should 
have pairwise distinct columns. Note that cases n = j (mod 3), j = 0, 1, 2, 
are considered by the same way. Suppose, say, n = 3t. According to (19.31) . we 
are interested in only cases when after removing r > 1 columns of A, all row 
sums will be odd. Suppose that after removing r columns of A, we have that 
p sums remain to equal to 3 and n — p sums equal to 1. This means that the 
total number of the removed l's equals to 2(n—p) = 6t — 2p. Since, removing 
a column, we remove three l's, then the number of the removed columns 
equals to r = 2t — y . Thus p = 3m and r = 2(t — m), m — 0, 1, t — 1. 
However, if m — t — 1, then r = 2. By the condition, these two columns 
are distinct, therefore, we conclude that at least one row sum equals to 2. 
The contradiction shows that the testing sequence is r = 4, 6, 2t. In cases 
n = 1, 2 (mod 3) we obtain the same testing sequence. 

Example 7. Let us check the parities of values of the permanent of circu- 
lants in Ay C A 7 . 

In this case t = |_f J = 2 and, therefore, the testing sequence contains 
only term r = 4. Note that matrix A r has all odd rows if and only if one 
row sum equals to 3 and each of 6 other row sums equals to 1. Indeed, let 
after the removing 4 columns of A, remain p sums equal to 3 and 7 — p sums 
equal to 1. This means that the total number of the removed l's equals to 
2(7 — p) and the number of the removed columns equals to r = 4 = 14 ~ 2p , 
i.e., p = 1. Moreover, since in a circulant all rows are congruent shifts of 
the first one, it is sufficient to consider the case when precisely the first row 
sum equals to 3 and others equal to 1 (the multiplication on 7 does not 
change the parity of the result). This opens a possibility of a momentary 
handy test on the parity every circulant of class Aj. This test consists of 
the removing all four columns beginning with 0. If now every rows 2, 7 
has one 1, then the permanent is even; otherwise, it is odd. We check now 
directly that from (^) = 35 circulants exactly 21 ones have odd permanent. 

Remark 1. In 1967, Ryser [13] did a conjecture that the number of the 
transversals of a latin square from elements l,...,n ( i.e., the number of 
subsets of its n pairwise distinct elements, none in the same row or column) 
has the same parity as n. If n is even, then the conjecture has been proved 
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by Balasubramanian [lj . Besides, in pQ Balasubramanian did a conjecture 
for the parity of a sum of permanents, such that the truth of this conjecture 
yields Ryser's hypothesis for odd n. In the same year (1990), using the result 
of Example 7, the author disproved Balasubramanian' s conjecture (private 
communication to Brualdi). It is interesting that soon Parker (see[5], p.258j 
indeed found several latin squares of order 7 with even number of transver- 
sals. Add that later (\18\) we found even an infinite set of counterexamples 
to the Balasubramanian conjecture. 

10. Open problems 

1. To prove the MCI (Section 7). 

2. (Cf.pH, PP- 171-172). Consider class A n (l, 1 + a, 1 + b), where < a < 
b < ^jp^. Since A n (l,l,l) = A 3 , then Voorhoeve's lower estimate for the 
permanent (12.161) trivially holds for matrices in A„(l, 1 + a, 1 + b). It is clear 
that, for a > 0, b > 0, it should exist an essentially stronger lower esti- 
mate. However, using Van der Waerden-Egorychev-Falikman theorem to 
class A ra ( 3+ ^ +fe , 3 b , 3 _^ b ) of doubly stochastic matrices, for the perma- 
nent of A n (l, 1 + a, 1 + 6)-matrices we obtain even weaker lower estimate of 
the order Ciy / n( 3+ " +b ) n << C(|) n . The problem is to find a stronger lower 
estimate for the permanent in A n (l, 1 + a, 1 + b). 

3. (Cf. [17], pp. 115-116). Let M be a circulant of order n with integer 
elements. We conjecture that, for every integer m, we have perM = 
(—l) n per(mJ n — M) (mod n), where n x n-matrix J n consists of l's only. 
A special case of this conjecture, for m = 1, M = I n + P + ... + P k ^ 1 
in the equivalent terms was formulated by Yamamoto [27] and proved for 
k < 3. The author [15] proved the truth of the conjecture in case m = 1 
for arbitrary circulant M ( including Yamamoto's conjecture for every k). 
In [T7J the conjecture was proved for every m and prime n. The question is 
open in case of composite n even in case k = 3. 

4. Two Latin rectangles let us call equivalent, if the sets of their elements 
in the corresponding columns are the same. Note that numbers |A^| one 
can treat as the numbers of equivalence classes of Latin triangles. Let 
A = I n + P + P 2 . In [20] the author proved that the cardinality of the 
corresponding equivalent class is 2 n + 6 + 2(— l) n . To find the cardinality of 
the equivalent class which is defined by matrix I n + P + P 3 . 
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